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Abstract 



We developed a 12-liter volume neutron detector filled with the liquid scintillator 
EJ301 that measures neutrons in an underground laboratory where dark matter and 



& neutrino experiments are located. The detector target is a cylindrical volume coated 

• iH on the inside with reflective paint (95% reflectivity) that significantly increases 

the detector's light collection. We demonstrate several calibration techniques using 

_C point sources and cosmic-ray muons for energies up to 20 MeV for this large liquid 



scintillation detector. Neutron-gamma separation using pulse shape discrimination 
with a few MeV neutrons to hundreds of MeV neutrons is shown for the first time 
using a large liquid scintillator. 
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1 Introduction 



The Sanford Underground Research Facility (SURF) was chosen as a site for 
ultra-low background experiments. The current two experiments are the di- 
rect detection of dark matter utilizing xenon with the Large Underground 
Xenon (LUX) [TJ experiment and the search for neutrinoless double-beta us- 
ing germanium with the Majorana Demonstrator [2] . Understanding the 
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in-situ background levels is extremely important for these rare-event physics 
experiments. Although positioning experiments in a deep underground labo- 
ratory significantly suppresses the background caused by cosmic-ray muons, 
the residual muons still create fast neutrons [3]. The intensity of the muon- 
induced neutrons depends largely on the depth of the underground labora- 
tory |3J. The energy spectrum, multiplicity, and angular distribution of the 
muon-induced neutrons are not particularly well measured. In addition, there 
are also fast neutrons from (a,n) reactions that are created in the surrounding 
rock by natural radioactivity. In this paper we demonstrate a neutron back- 
ground characterization technique that has been developed for a large liquid 
scintillation detector. 

One major problem in neutron detection is the separation of neutrons from 
the electromagnetic background caused by gamma rays from the environment 
and internal contamination of the detector materials. The pulse shape discrim- 
ination technique, which uses the difference in the shape of the scintillation 
pulses generated by neutrons and gamma rays, has been implemented success- 
fully with small neutron detectors for many years [If5|6|7f8] . Unfortunately, 
small neutron detectors (of a few liters or less) have low efficiency when de- 
tecting neutrons and are high cost compared to size and efficiency For the 
successful neutron background characterization, we require a higher efficiency 
of neutron detection and a broader energy sensitivity, along with lower cost 
per detector. The development of a relatively large volume detector increases 
the possibility of significantly improving the neutron detection efficiency. More 
importantly, it opens a window for exploring the neutron energy at a few MeV 
up to a few hundred MeV. Although most of the neutron-gamma discrimina- 
tion experiments were carried out using small detectors, the possibility of 
neutron-gamma discrimination using time of flight (TOF) measurements have 
also been investigated using a large neutron detector [9] . In this work we inves- 
tigate the possibility of using a large volume detector for the direct detection 
of neutrons with energy of a few ten MeV, with pulse shape discrimination to 
distinguish from gamma rays. 



2 Experimental setup 



We have constructed a liquid scintillation detector that is 1 meter long and 
5 inches (12.7 cm) in diameter. This detector is fabricated using a aluminum 
cylindrical housing with two Pyrex windows on each side of the cylinder at- 
tached to PMTs. To mitigate light loss, the inner surface of the detector was 
covered with reflective paint EJ520 (Eljen Technology) that has 95% reflectiv- 
ity. The reflective paint makes scintillation photons scatter multiple times from 
the detector walls, which partially compensates for the relatively poor PMT 
photo-cathode coverage. The detector volume was filled with liquid scintillator 




Fig. 1. The setup of our liquid scintillation detector. The aluminum foil pans are 
required for safety in case of leaks. 

EJ301, which is specifically designed for neutron-gamma discrimination. Scin- 
tillation light is collected by two 5-inch Hamamatsu R4144 PMTs attached 
to both Pyrex windows. The geometry of the detector setup is presented in 
Fig. [I] The total photo-cathode coverage of the detector is at about 3% and 
the PMT quantum efficiency is about 20% for photons with a wavelength of 
300-700 nm. The EJ301 is composed of carbon and hydrogen atoms, and has a 
H/C ratio of 1.212 and a density of 0.874 g/cm 3 . Its light output is 78% as of 
anthracene with a maximum emission at 425 nm, which exactly corresponds 
to the most sensitive region of the Hamamatsu PMTs. 

Contamination of oxygen in the liquid scintillator results in a reduction of 
light output. In order to minimize the amount of oxygen contamination, the 
scintillator was thoroughly purged with dry argon and then the detector vol- 
ume was sealed before the measurements were performed. The operational 
voltages of both PMTs were determined to be 2000 V. In order to avoid satu- 
ration from high energy events, especially those close to the PMTs, both PMT 
output signals are attenuated at 23 dB before connecting to the DAQ. The 
DAQ consists of a fast flash ADC that analyzes the PMT output signal. The 
sampling frequency of the flash ADC is 170 MHz, providing a signal amplitude 
every 5.88 ns. The ADC is controlled by a program based on the MIDAS data 
acquisition system software [TO] . 



3 Energy and position calibration procedures 



For a large scintillator, the detector response to energy with each individ- 
ual PMT is determined by the energy deposition and the position where the 
incident particles interact with a target in the scintillator. This position de- 
pendence is caused by the attenuation of the light in the tube [UJ. To obtain 
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Fig. 2. Position of the incident particle and the total light yield being split and 
collected by individual PMTs. 

the position information, the charge ratio from the two PMTs is used to char- 
acterize the position of an incident particle. The following criteria are applied 
to select events: 

(1) Both PMTs must be triggered. 

(2) Both PMTs must not be saturated. 

(3) Time coincidence is within 30 ns, which is the time difference of the 
largest sample in the pulse between two PMTs. 

Assume an energy deposition of E tot is created at a distance of X (distance to 
the middle, see Fig. |2|), the light collection by two PMTs Li e f t and L rig h t , and 
the total light yield L tot proportional to E to t- Without any energy loss, the total 
light yield would be evenly split by the two PMTs, L; e / t = L rig ht = 0.5L tot . In 
reality we have to consider light attenuation, especially for such a big detector. 
Taking I as the attenuation length in the scintillator and D as the total length 
of the tube, a simple calculation in Eq.(fl| shows that the position can be 
determined by the combination of L/ e /j and L rig ht, respectively. 




■■ 0.5L tot e 
0.5L tot e 

■ X/l. 



-{D/2-X)/l 
-(D/2+X)/l 



(1) 



The light collection Li e f t and L r i g ht will then be converted to photoelectrons 
at the photocathode of the PMTs. If we define a0 to stand for the total 
charge converted from Lj e /t and ol the total charge converted from L r i g ht, 
aO should be proportional to L te f t (a0 oc Li e f t ) and al proportional to L right 

(ol oc L r i g ht)- Thus, we can use the charge ratio, In JaO/al oc X/l, to interpret 
the position of the particle along the tube. 

The detector response to energy with respect to one of the PMTs for a surface 
background run is presented in Fig. |3J The value of X/l varies from -2 to 
2 which is four times the attenuation length from one end of the tube to 
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Fig. 3. Detector response to energy with respect to one of the PMTs for a sur- 
face background run. X/l stands for the relative position of the particle. aO is the 
integrated charge (the pulse area) of the PMT0. 

the other. Considering the length of the tube is 1 meter, a mean attenuation 
length of ~25 cm is determined for the light transport in the scintillator. Note 
that this result is much less than the expected attenuation length for a EJ301 
scintillator (I > 1 meter). This is because the the light transport in the tube 
is dominated by the diffusive reflection, which dramatically exacerbates the 
attenuation. 

In the lower energy range, the detected events are dominated by gamma 
rays from the internal contamination of the detector components. In contrast, 
events at the high energy range are mainly from the cosmic muons. Clearly, 
the higher energy curve in Fig. [3] indicates a muon minimum ionization peak of 
~20 MeV at a 5-inch diameter. This minimum ionization has also been verified 
using a GEANT4 simulation [12]. The incident cosmic ray muons should be 
uniformly distributed along the tube (except for edge effects at both ends of 
the tube), which provides us a natural energy calibration source for energy up 
to 20 MeV. This curve shows that the energy response (total charge collected 
by PMT) of a single PMT has strong position dependence. Within the tube, 
i.e., X/l G [—2,2], there are three blocks filled with no data. The blank area 
at the bottom in Fig. [3] is caused by the energy threshold set for the trigger 
of the PMTs while the blank areas on both sides are caused by the saturation 
of the PMTs. 
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Fig. 4. The position (left) and energy (right) calibration from a 22 Na source in 
terms of the 1.275 MeV gamma ray line. The data in the left plot is fitted by a 
third-order polynomial function. 



Energy-scale calibration is performed using a 22 Na gamma ray source. The 
22 Na radiation source produces two gamma ray lines with 0.511 MeV and 
1.275 MeV energies. Measurement is performed by setting the uncollimated 
source on the top of the aluminum tube every 2.5 cm from one end to the other. 
Calibration is made using the 1.275 MeV gamma ray line from a 22 Na source. 
The left plot of Fig. [I] gives a position calibration interpreted by X/l. It results 
in nearly linear relation between X/l and the actual position. The right plot 
shows a position-dependent energy response to the 1.275 MeV gamma rays 
from the source. 



The position dependence of energy makes the energy calibration complicated. 
Since both ends of the aluminum tube are coupled to PMTs, the energy scale 
can be position independent [13] based on the Eq.Q. 
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Since J Li e f t L rig ht is proportional to y/a0 x al, we can use \/a0 x ol to rep- 
resent the energy scale and avoid the position dependence along the tube. To 
verify this hypothesis, a plot of y/aO x al versus X/l has been created in Fig. 
[5] where a distinct and almost horizontal line is shown for the muon minimum 
ionization process with the y/a0 x al at ~ 4000. Other than the edge effect at 
both sides, it shows that the position independence of energy is quite accurate 
if we use y/a0 x al to describe the total energy deposition. Using 22 Na and 
AmBe [H] sources we also obtain consistent results for gamma ray lines at 
1.275 MeV and 4.4 MeV. Therefore, an energy calibration curve is obtained 
for few MeV up to 20 MeV. 
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Fig. 5. Plot of y/aO x al vs X/l to show the position independence of energy along 
the tube. This is a surface run where an AmBe source was placed on top of the tube 
and 22.5 cm from the right end. 

4 Neutron-gamma discrimination 



Measurements of neutron-gamma discrimination were carried out using an 
AmBe neutron source [13] • It was placed vertically 6.5 cm directly above the 
outer wall of the tube and horizontally 22.5 cm away from PMTO. The neutron 
emission rate of this source is ~100 Hz, and the neutron energy spectrum is up 
to 11.2 MeV. The digitized pulses from each event have been recorded for data 
analysis. The analysis procedure includes calculating the total charge per pulse 
per event and the charge that corresponds to the tail of the pulse (delayed 
charge). The total charge is defined as the integral under the pulse from 8 
samples before to 40 samples after the peak, where the pulse sample is defined 
as the signal amplitude every 5.88 ns. The delayed charge is the integral under 
the pulse from 8 to 40 samples after the peak. Since the decay-time constants 
for scintillation light from neutrons and gamma rays are different, the neutron- 
gamma discrimination can be obtained using the ratio of the delayed charge 
to the total charge, -4^y, where d(A) is the delayed charge integrated over the 
delayed pulse area and t(A) is the total charge integrated over the total pulse 
area. 



A dedicated study shows that the neutron-gamma separation can only be 
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Fig. 6. The performance of neutron-gamma discrimination from PMTO. Relative 
position X/l at 0, 0.5, 1 and 1.5 are displayed in the individual plots. 

found in a very narrow range of position along the tube (see Fig. pj). 

The plots describe the separation of neutrons and gamma rays in terms of 
their pulse shape differences. It clearly shows that neutron-gamma separation 
decreases when energy deposition occurs far from the target PMT. This is 
expected because the diffusive reflection, absorption and re-emission of light 
washes out the pulse shape difference between neutron and gamma ray events, 
which happen farther away from the target PMT. Therefore the charge ratio 
delayedl/al from PMT1 is adopted to separate the neutrons at one end of 
the tube (X/l < 0) and delayedO/aO from PMTO is used for the other end 
(X/l > 0). The neutron-gamma separation vanishes for a broad position range 
as shown in Fig. [71 

The plots in Fig. [6] indicate that the gamma ray band and the neutron band are 
separated when position is held constant. However, they have different slopes 
as the position varies, which results in the separation disappearing when we 
use a wide position range. Since y/aO x al represents energy, which is found 
to be independent of the position, the energies (x-axis) of the invariant data 
points are kept and the value of y, or ratios of the charge between the delayed 
and the total are projected to the horizontal line for all data points (according 
to the slope k, i.e., (x, y) =>■ (x,y — kx)). The relation between the position 
and the slope of the gamma ray band is interpreted in Fig. [8j 
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Fig. 7. Neutron-gamma discrimination for a broad position range. The top two 
plots adopt the signals from PMT0 at position ranges X/l > and X/l > 0.5 while 
the bottom two plots take the signals from PMT1 at position ranges X/l < and 
X/l < -0.5. 

After summing all the narrow position slices together, a combined neutron- 
gamma separation plot is obtained as shown in Fig. [9j In order to get better 
neutron-gamma separation, the condition \X/l\ > 0.5 is required to limit 
the positions to values that are not too far away from the PMTs. In order 
to compare the improvement of the separation, the figure of merit [IS] that 
describes the goodness of separation between gamma rays and neutrons is 
defined below: 



Goodness 






(3) 



where R is a mean of the distribution for the ratio between the delayed charge 
and the total charge and a is the Half Width at Half Maximum (HWHM) of 
the distribution of the charge ratio. The left plot in Fig. [10] illustrates how the 
goodness of separation is calculated for the position range X/l < —0.5 and 
the energy at 40 MeV. A Gaussian distribution is assumed for both neutrons 
and gamma rays. The multiple peaks in the neutron band indicate different 
recoils induced by protons, alphas, and 12 C The right plot in Fig. 10 gives the 



comparison of the figure of merit before and after the position cut, which is 
illustrated in Fig. [9j It shows that the neutron-gamma separation gets better 



Q. 
O 

C/) 



10" 4 - 



10 5 - 




-1.5 



-1 



-0.5 



0.5 



1 



1.5 



2 
X/l 



Fig. 8. The relation between the position and slope of the gamma ray band, 
after position cuts are applied. 

Given the neutron energy spectrum of the AmBe source, the maximum neu- 
tron recoil energy detected by the detector is up to 6 MeV. However, for a 
surface run, the cosmogenic neutrons can create higher energy events. Because 
the AmBe source is much closer to PMTO than PMT1, we see more source 
related events in the upper two plots of Fig. [9] There are distinct peaks in 
the gamma-ray band, which are the muon minimum ionization peaks for ~ 20 
MeV and 4.4 MeV gamma-ray lines from AmBe source. Although the muon 
ionization peak looks diffused in a wide range due to the edge effect, it gives 
a well defined peak once we remove the range that is close to the PMT ends. 
The neutron-gamma separation performs well up to 100 MeV, which indicates 
this detector is able to detect neutrons with energies of up to a hundred MeV. 
The low energy range (< 6 MeV) of its neutron band is overwhelmed by AmBe 
neutrons. However, this cannot be separated from the gamma-ray band due 
to the contamination of 4.4 MeV gamma rays from the AmBe source. This is 
because some gamma rays associated with neutrons enter into the detector at 
the same time. In order to remove this source related gamma-ray contamina- 
tion, 4 inch lead bricks were added between the AmBe source and the detector. 
The whole setup is simulated using GEANT4 and its geometry is displayed in 



Fig. 11 



Rather than implementation on the surface, the new AmBe run was performed 
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Fig. 9. The corrected neutron-gamma separation results in terms of signals from 
PMTO (top two) and PMT1 (bottom two). The energies are also corrected based 
on calibration. 
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Fig. 10. The left plot is a projection on the charge ratio, delayedl/al, for the 
position range X/l < —0.5 and the energy at 40 MeV. The right plot shows the 
comparison of the Figure of Merit that characterizes the goodness of separation 
between neutrons and gamma rays versus the visible energy. 

underground at Soudan mine (2100 m.w.e). This has the extra benefit of dra- 
matically reducing the high energy backgrounds from cosmic rays. The entire 



process of data analysis and reconstruction is shown in Fig. 12 Comparing 
the upper right plot in Fig. [9] to the upper left plot of Fig. [12J one can clearly 
see that the high energy component has almost disappeared. This is expected 
because the 800 meters of rock overburden highly suppresses the intensity of 
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Fig. 11. Setup for the run with an AmBe source at Soudan Mine underground. The 
source is located 22.5 cm away from the right hand PMT and above two layers of 
lead bricks (total 4 inch). 
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Fig. 12. A calibration with an AmBe source at the Soudan mine underground. 
The upper left plot is a combination of the AmBe source signal and background 
data. The upper right plot gives source signal with background subtracted. The 
bottom left plot shows neutrons separated from gammas at 1.3 MeV energy and for 
a delayedO/aO < 0.25 charge ratio cut. The bottom right plot displays the processes 
for neutron energy reconstruction. The result is compared with simulation. 
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cosmic rays from the atmosphere. After the background is subtracted, as seen 
in the upper-right plot of Fig. [12J the low energy gamma rays from the in- 
ternal contamination are thoroughly removed. The background subtraction is 
not very effective at the high energy muon tail due to the limitation of our 
statistics. The 4.4 MeV gamma ray line from the AmBe source is detected 



by the detector shown in Fig. [12} A 1.3 MeV energy cut is set on the neu- 
tron band to cut off data that is affected by random noise at low energy. 
This 1.3 MeV visible energy is equivalent to 4 MeV recoil energy by applying 
the quenching factor measured for EJ301 (NE213) and BC-501A by several 
groups [T6fl7|ll8] including KamLAND |19j . Quenching factor is a convolution 
of the ionization efficiency and scintillation efficiency induced by nuclear re- 
coils [20J. These two efficiencies are well described by Lindhard's theory [21] 
and Birks law |22j. Although the composition of the liquid scintillator we 
use is a little bit different from that of KamLAND, both quenching factors 
induced by nuclear recoils are dominated by nearly the same ionization effi- 
ciency. Therefore, a small difference of less than 5% in the quenching factor 
between EJ301 (NE213) and KamLAND's scintillator offers a cross check and 
confirms our use of the quenching factor. Pure neutrons are obtained after 
appropriate separation from gamma ray and noise as shown in the bottom left 



plot of Fig. 12 In order to determine the recoil energy, the quenching factor 



is applied to the visible energy of the neutrons, as shown in the bottom right 



plot Fig. 12 A GEANT4 based simulation has been performed to compare 
with the measurement data. A good agreement is found for recoil energies 
above 4 MeV. 

The neutrons and gamma rays are not well separated around 4 MeV, as shown 



in the upper right plot of Fig. [12} This indicates that 4 inch thick lead bricks 
might not be thick enough to block such high energy gamma rays. The neutron 
signals are still suffering contamination from gamma rays, although there is 
an improvement over no shielding. In addition, the low energy range of the 
neutron band has an overlap with the gamma ray band. This is because the 
pulse shape of low energy events is easily affected by the attenuation of light 
in the tube and the fluctuation of the pedestal level of the ADC. A 23 dB 
attenuator applied to the signal makes the signal-to-noise ratio even worse. The 
trigger rate is very high in the low energy range, which was mainly contributed 
as internal contamination in the detector. It will become a challenge if we 
attempt to measure low energy (a,n) neutrons in an underground laboratory. 
We may have to consider sacrificing the sensitivity at low energies to improve 
the data quality. 

A test has been conducted using the same AmBe run, but with a raised the 



energy threshold as shown in Fig. 13 (maxSam = maximum sample - pedestal) 
for both PMTs. The top two plots show those events with low energies and 
the positions close to one PMT. This causes a weak signal at the other PMT 
that is removed by raising the threshold. These events are poor quality in 



13 




Visible Energy (MeV) 



Visible Energy (MeV) 




1 



O 0.4 
Tj 0.35 - 

cd : 

I "E- 

■D 
0.25 

0.2 E- 

0.15 
0.1 
0.05 F- 



20 25 30 

Visible Energy (MeV) 








15 20 25 30 

Visible Energy (MeV) 



Fig. 13. The improvement of data quality from raising the energy threshold 
for both PMTs. The left three plots are the data collection with the conditions: 
maxSamO > 30, maxSaml > 30; the right three plots are the results with cuts: 
maxSamO > 110, maxSaml > 130. The bottom four plots apply a further cut on 
X/l > 0.5. 

terms of their pulse shapes to at least one of the PMTs. Since the value of 
maxSam represents the magnitude of the energy deposition in the same way 
as the integrated charge, the ratio of (a0/al) : (max SamO/ maxSaml) should 
be a constant if both PMTs are triggered by the same event. 



The middle two plots in Fig. [13] clearly show how low energies are affected by 
poor quality events and how they have been improved by raising the energy 
threshold. The ratio, (aO/al) : (max SamO/ maxSaml) is found to be a con- 
stant (~1) from 2 MeV up to a few ten MeV. The comparison of the bottom 



two plots in Fig. 13 gives more information. First, the phenomenon of neutron- 
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gamma superposition in the low energy range is eliminated. Second, neutrons 
can be easily identified from the gamma ray band. Finally, gamma rays are 
well characterized by themselves, e.g. the 2.6 MeV gamma ray line from ra- 
dioactive decay and the 4.4 MeV gamma-ray line from the AmBe source. The 
high energy events are barely affected by raising the threshold. It is worth 
pointing out that the energy threshold of the neutron band is slightly larger 



than the gamma ray band in the bottom right plot Fig. 13 This is because we 
set the same threshold on the maximum sample for both neutrons and gamma 
rays. However, the delayed pulses induced by neutrons is larger than that of 
gamma rays. Therefore, the energy of neutrons represented by the integrated 
area is slightly larger than that of the gamma rays. 



5 Discussion and Conclusion 



Other than protons, high energy neutrons induce recoil deuteron, alphas, and 
12 C in scintillators, which causes distinct bands in the pulse shape spectrum for 
small size detectors [23]. For a large detector, these features are also observed 
at positions very close to the PMTs. For those recoils occurring away from 
PMTs, the features are washed out due to the complexity of light transmission. 

With a 12 liter liquid scintillation detector, the surface muon minimum ion- 
ization peak is utilized to calibrate the detector with an energy of up to ~20 
MeV. The representation of position using the parameter X/l is found to be 
very convenient for interpreting the features of light transport in the scintil- 
lator. Energy independence, in terms of position, is found based on a simple 
mechanism, which makes the energy reconstruction much easier. It is also 
useful to separate neutrons from gamma rays when a wide range of positions 
are combined, because of this, we have demonstrated the position dependent 
neutron-gamma separation. A pulse shape analysis procedure has been per- 
formed to distinguish neutrons from gamma rays. A new algorithm has been 
developed for a large liquid scintillation detector to directly measure neutrons 
at a few MeV to a few hundred MeV. 

The direct measurement of fast neutrons, traditionally, can be performed using 
the time of flight (TOF) technique. However, the small solid angle coverage 
dramatically limits its efficiency Because of the extremely low intensity of 
neutrons, it is almost impossible to perform such a measurement for the fast 
neutrons in a deep underground laboratory. Rather than using TOF measure- 
ments, we provide a new method to directly measure high energy neutrons 
with a much better detection efficiency (~30% at 10 MeV for the AmBe run) 
and at a relatively low cost. Depending on the needs of the underground depth, 
an array of such neutron detector modules could be employed to increase the 
detection efficiency. We conclude that a large scintillation detector can be used 
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to measure fast neutrons for ultra-low background experiments underground. 
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